Statistical properties of one dimensional attractive Bose gas 
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Using classical field approximation we present the first study of statistical properties of one di- 
mensional Bose gas with attractive interaction. The canonical probability distribution is generated 
with the help of a Monte Carlo method. This way we obtain not only the depletion of the conden- 
sate with growing temperature but also its fluctuations. The most important is our discovery of 
a reduced coherence length, the phenomenon observed earlier only for the repulsive gas, known as 
quasicondensation. 
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INTRODUCTION 



Statistical properties of finite number of bosons con- 
fined by a trapping potential are intensively studied 
ever since the very first experiments on dilute gas Bose- 
Einstein condensates. Initially a full understanding has 
been reached of the statistics of the confined ideal Bose 
gas PHI)]. In this context criticism of the most widely 
used grand canonical ensemble has been raised. As it 
turns out it predicts unphysical fluctuations of the num- 
ber of condensed atoms. Interactions, an essential as- 
pect of Bose atoms physics make the problem of statistics 
rather complicated. In the simplest, although still aca- 
demic case, that of the box with periodic boundary condi- 
tion, at least we know a priori what is a condensate. It is 
the zero momentum component of the atomic field. In the 
box with periodic boundary conditions we know more: we 
know analytically the spectrum and the expressions for 
the collective Bogoliubov excitations. So, one can look at 
the statistics of the gas as that of noninteracting bosonic 
quasiparticles. This determines the statistical properties 
of the gas at least at low temperatures Q. Extension 
of the method to higher temperatures is already signif- 
icantly more complicated. The Bogoliubov-Popov spec- 
trum depends on the actual number of the condensed 
atoms rather than on their total number. Thus, itself is 
a fluctuating variable. Several papers attempted to take 
this into account H Q . Another serious problem is the 
assumption that quasiparticles do not interact. In fact 
they do and they have a finite (temperature dependent) 
lifetime. In a realistic harmonic trap we have no analytic 
formulae for the quasiparticles and also the condensate 
degree of freedom is known only a posteriori since it is 
a function of an interaction strength, a frequency of the 
harmonic trap, a number of atoms and a temperature of 



the sample. In two recent papers [10L we have shown 
how to overcome most of these problems. To this end we 
proposed to use the classical field description of the sys- 
tem generating the proper thermal equilibrium distribu- 
tion using a classical Metropolis algorithm 1121 . Details 
of classical field approaches are presented in [13, |l4| . It 
is the purpose of this Letter to report the first study of 
statistical properties of an attractive Bose gas. Confined 
to a box it is unstable. Hence, this simple model makes 
no sense for the attractive gas. The situation is different 
in a harmonic trap. Limited attractive condensates do 
exist in two and three dimensions. The lithium- 7 con- 
densate was among the first ones to be observed [lal . It 
then has lead to creation of bright solitons [3, [TtJ ■ In 
strictly ID systems the effective repulsion of the kinetic 
energy makes it stable for any number of particles. We 
therefore concentrate our attention on the one dimen- 
sional, finite, attractive Bose gas trapped in a harmonic 
potential. Thus the Hamiltonian of the system has the 
form: 
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where \fr is the atomic field operator, m is the mass of 
an atom, loq is the frequency of the harmonic trapping 
potential and we are particularly interested in g < 0. The 
case of positive coupling constant was studied in detail 
in [IH. 

The classical fields approximation consists in replace- 
ment of the quantized atomic field by a c-number wave 
function. This wave function is then conveniently ex- 
panded as a sum over the harmonic oscillator wave func- 
tions and expansion coefficients a n are classical stochas- 
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tic complex variables. 
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where the oscillator eigenfunctions y n are chosen to cor- 
respond to a harmonic oscillator of the frequency uj. Note 
that uj is not necessarily equal to ojo- 
Thus the energy functional: 

E[*\=j **(x)(^ + \mu%A*{x) + \ J \^{x)\\ 
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takes a form: 
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with the interaction energy £?j„t being a quartic form in 
the variables a n and (n|a; 2 |n'} is the matrix element in 
the oscillator basis. 

We are going to look at the statistics of a sample com- 
posed of N atoms, so the amplitudes are subject to a 
constraint: 



n=0 
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Having converted a quantum statistical physics problem 
into classical one we then define a canonical equilibrium 
distribution of the amplitudes a n as 



P oc exp 



(6) 



An efficient algorithm generating this probability distri- 
bution has been invented by Metropolis (l2^ . 

It is easy to check with the help of Gross-Pitaevskii 
equation imaginary time propagation that the zero tem- 
perature attractive ID Bose gas has a wave function 
which is very close to a Gaussian. Its width shrinks with 
growing number of condensed atoms and with growing 
absolute value of the negative coupling constant g. This 
is our guide to the optimal choice of the frequency to 
defining the actual oscillator base in the expansion of the 
atomic classical field © . Thus this frequency becomes a 
function of the coupling and of the temperature as they 
both determine the width of the condensate wave func- 
tion. More precisely: as we have shown in [13] the clas- 
sical fields reproduce faithfully the available analytically 
exact probability distribution of a harmonically confined 
ideal Bose gas for the cut-off parameter n max satisfying: 




FIG. 1: (Color online)Effective frequency uj(N,T) for various 
values of the parameter g. The total number of atoms N = 
500. Solid red line corresponds to g = —0.01, dashed green 
to g — —0.0067 and dotted black to g = —0.003. In inset we 
compare the FWHM of the Gaussian with effective uj with 
the width of the condensate (crosses). 



where T is the temperature of the gas. A natural gen- 
eralization of this condition for a weakly attractive ID 
Bose gas is 



c frw{N,T) = k B T 



(8) 



where the frequency u>(N,T) is chosen such that the 
corresponding Gaussian function describes the ground 
state of the Gross-Pitaevskii equation for No interacting 
atoms. Strictly speaking No should be chosen self consis- 
tently. In this Letter, however, we are satisfied with No 
being a number of condensed ideal gas atoms at given 
temperature. The general rule for determination of the 
condensate wave function follows from the Onsager- Pen- 
rose pj| definition calling for the diagonalization of the 
one-particle density matrix 



Pi,j = (a*Uj) = ^KPt{n)Pj(n) 



(9) 



,hoj = k B T 
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and identification of the condensate wave function with 
the eigenvector corresponding to the leading eigenvalue. 
As a cross-check we verified that our identification of base 
indeed yields the condensate wave function very close to 
the n=0 state. 

It is clear that our choice of the frequency oj(N,T) 
makes it monotonically decrease with the temperature 
and tend to that of the empty trap: loo- In Fig. Q] we 
illustrate this for several values of the coupling g. In the 
inset we also confirm a consistency of our choice of the 
base by comparing the width of the ground state of the 
Gross-Pitaevskii equation with the width of the actual 
condensate wave function computed via diagonalization 
of the one- particle density matrix. Small steps visible 
in the inset result from taking the solution of © for the 
cut-off as the nearest integer. 
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FIG. 2: (Color online) Number of atoms in the condensate as 
a function of the interaction strength g for two fixed tempera- 
tures T = 0.26Tb (squares) and T = 0.63Tb (crosses). These 
temperatures are indicated in the inset, where the decay of 
condensed atoms in the ideal gas is presented. 



FIG. 3: (Color online) Dispersion of number of condensed 
atoms as a function of the interaction strength g for two 
fixed temperatures T = 0.26Tc (squares) and T = 0.63Tc 
(crosses). These temperatures are indicated in the inset, 
where the fluctuation of condensed atoms in the ideal gas 
is presented. 



Multiparticle Hamiltonian with interactions dependent 
on the mutual distances between atoms with the har- 
monic potential trapping has the center of mass of the 
system obey interaction free Schrodinger equation de- 
coupled from all other degrees of freedom. Thus this 
variable is entirely missing from any mean field approx- 
imation in which a nonlinear Gross-Pitaevskii equation 
describes the system. This does not have important con- 
sequences for the repulsive gas. In this case, the conden- 
sate is broader than the ground state of the harmonic 
potential, so the uncertainty over the position of the cen- 
ter is just a correction. The situation for attractive gas is 
different. Now, the condensate is spatially squeezed and 
the shot-to-shot variation of the position of the center of 
mass is at least as broad as the ground state of the empty 
harmonic potential 3, 2(3] • We therefore stress that our 
approach disregards this uncertainty, so the statistical 
properties derived here pertain to the reference frame of 
the center of mass of the system. Measurements in the 
center of mass system can meet some difficulties due to 
uncertainty of position of the center of mass. However, 
a single shot simultaneous detection of many particles 
resolves this problem [5(3] . 

In our classical field approximation effects of quantum 
fluctuations and quantum depletion are missing. The 
method is suitable for weakly interacting Bose gas in the 
quantum degenerate regime. Throughout this paper we 
use the oscillator units of position, energy and tempera- 
ture, y-^j-t fo^o an d ip 1 respectively. Hence a dimen- 
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sionless coupling g is in units of 

All presented results are for N = 500 atoms and the 
temperature is in units of the characteristic temperature 
indicating transition to the quantum degenerate regime 



of a finite sample of an ideal Bose gas N = Tcln(2Tc) 
where N is the total number of atoms 21). In Fig. [2] 



we plotted a number of condensed atoms as a function 
of the coupling constant g for negative and also positive 
values (using the method explained in detail in [11|) for 
two temperatures, one very low (T/Tq = 0.26) and the 
other T/Tq = 0.63. Since the condensation depends on 
the local density of particles, the effects of positive and 
negative interactions are opposite. The first one swells 
the condensate reducing the local density with respect to 
the ideal gas while the other shrinks it enhancing the con- 
densation process. Thus the number of condensed atoms 
is a monotonically decreasing function of the interaction 
parameter g as we go from the negative to the positive 
values. In Fig. [3] we present the relative fluctuations 
also as a function of changing sign g and for the same 
temperatures as in Fig. [2l General observation is that 
fluctuations of the condensed atoms number grow with 
the strength of the interaction in both negative and pos- 
itive directions [2tI | . An unexpected feature, however is 
that the minimal fluctuations are observed not for g = 
but for tiny attractive interaction of —0.0025. 

Perhaps the most interesting aspect of a ID repulsive 
Bose gas predicted [22j and then observed experimentally 



23. 



24) is the phenomenon of a quasicondensation. It is a 
reduction of the coherence length to values smaller than 
the length of the condensate above a certain characteris- 
tic temperature. In the recent Letter [SBJ the correlation 
length has been related to the gray solitons formed dur- 
ing a rapid cooling process [2^. We have computed the 
temperature dependent correlation length for the attrac- 
tive case and found that the quasicondensation occurs 
also in this case. This is probably the most important 
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FIG. 4: (Color online) The width of the first order correla- 
tion function and the width of the condensate as functions of 
temperature. 



result of this Letter. We note that one does not expect 
gray solitons to appear in rapidly cooled attractive Bose 
gas, so above mentioned connection to solitons certainly 
does not hold in this case. The first order correlation 
function is defined as 



9l(x,-x) 



(tf(-a;)**(a;)) 
<l*(z)| 2 > 



(10) 



In Fig. [5] we plot its full-width at half maximum as a 
function of temperature and compare it with the width 
of the condensate. The effect of quasicondensation is 
visible. It occurs at the temperature T m O.QTc- 

Summarizing: We have presented the first study of 
equilibrium thermodynamics of the attractive ID Bose 
gas trapped in a harmonic potential. Our results are for 
the canonical statistical ensemble, thus the temperature 
is a control parameter. The classical field approxima- 
tion is used. The appropriate probability distribution 
is obtained numerically using a Monte Carlo technique. 
We found a reduced coherence length above some char- 
acteristic temperature. This phenomenon was previously 
known only for a repulsive gas. 
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